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Abstract — Plasma  convection  at  high  latitudes  is  strongly  controlled  by  the  direction  of  the  interplanetary 
magnetic  field  (IMF).  When  the  IMF  is  southward,  two-cell  convection  patterns  exist  with  anti-sunward 
flow  in  the  polar  cap  and  return  flow  equatorward  in  the  auroral  oval.  When  the  IMF  is  northward, 
however,  the  situation  is  more  complicated  and  several  convection  patterns  have  been  proposed  to  explain 
the  observations  of  sunward  polar  cap  flow,  including  three-cell,  four-cell,  and  distorted  two-cell  patterns. 
Since  the  different  convection  patterns  are  expected  to  display  different  ionospheric  signatures,  the  Utah 
State  University  Ionospheric  Model  can  provide  insight  as  to  which  convection  model  applies  to  given 
situations.  The  model  was  run  to  obtain  diurnally  reproducible  lonosphenc  densities  and  temperatures  for 
summer  and  winter  conditions  using  both  distorted  two-cell  and  three-cell  convection  patterns.  Differences 
due  to  the  different  convection  patterns  manifest  themselves  in  the  depth  and  location  of  polar  holes  in  the 
F-region  electron  density.  While  the  total  depth  of  the  model  holes  is  a  characteristic  of  the  diurnally 
reproducible  pattern,  the  features  appear  and  are  recognizable  within  0.5  h.  Langmuir  probe  data  from  41 
DE- 2  passes,  during  which  the  IMF  B.  component  was  northward,  have  been  qualitatively  checked  against 
the  model  predictions.  The  cross  polar  cap  electron  density  profiles  of  a  large  majority  of  the  passes  more 
closely  conform  to  the  distorted  two-cell  runs  for  both  polarities  of  the  IMF  B,  component.  This  test  can 
be  generalized  to  rule  out  proposed  convection  patterns  based  on  the  presence/absence  and  position  of 
polar  electron  density  holes. 


INTRODUCTION 

Under  northward  interplanetary  magnetic  field  (IMF) 
conditions  electric  fields  in  the  polar  cap  (i.e.  the 
region  poleward  of  the  auroral  belt,  that  is  commonly 
characterized  under  southward  IMF  conditions  by 
dawn-dusk  electric  fields  and  anti-sunward  con¬ 
vection)  are  often  highly  structured  with  a  variety  of 
irregularity  dimensions.  The  irregularities  often 
include  structures  in  which  the  sign  of  the  dawn- 
dusk  electric  field  component  is  reversed,  producing 
sunward  convective  motion.  Large-scale,  distinct 
regions  in  which  the  fields  are  predominantly  of  one 
sign,  and  which  may  be  identified  with  convection 
cells  may  or  may  not  be  present. 

Previous  models  of  convective  electric  fields  for 
northward  IMF  have  proposed  patterns  with  variable 
numbers  of  cells.  Crooker  (1979),  in  studies  of  allowed 
merging  regions,  postulated  a  large  convection  cell 
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associated  with  the  polar  cap  and  the  lobes  of  the 
magnetosphere  for  IMF  B.  North  conditions. 
Circulation  in  the  cell  changed  direction  with  IMF  i 
For  strong  northward  IMF,  the  cell  was  concentrated 
on  one  side  or  the  other  of  the  polar  cap  depending 
on  the  size  of  Br  The  concept  of  a  four-cell  pattern 
was  introduced  by  Burke  et  al.  (1979)  to  explain  the 
53-2  satellite  observation  of  sunward  convection  in 
the  center  of  the  polar  cap.  Reiff  and  Burch  (1985) 
proposed  a  four-cell  pattern  with  two  lobe  cells  in  the 
polar  cap  and  two  narrow  cells  at  adjacent  lower 
latitudes  driven  by  viscous  interaction  in  the  boundary 
layer.  One  lobe  cell  could  dominate  the  other  depend¬ 
ing  on  the  strength  and  polarity  of  the  By  component. 
For  strongly  positive  B„  one  of  the  lobe  cells  in  each 
cap  was  driven  while  the  other  was  a  reclosure  cell  on 
closed  field  lines  driven  by  merging  in  the  opposite 
polar  cap.  Heelis  et  al.  (1986)  also  postulated  four¬ 
cell  patterns  with  the  dominant  lobe  cell  being  depen¬ 
dent  on  the  By  component.  Banks  et  al.  (1984)  and 
Clauer  and  Banks  (1936)  developed  convection 
patterns  having  four  cells  that  were  related  to  the 
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dayside  cusp  region  based  on  direct  coupling  of  solar 
wind  electric  fields  into  DPY  and  DPZ  related  current 
systems.  Their  model  resulted  in  two  extra  small  cells 
appearing  in  the  vicinity  of  the  cusp,  but  not  extending 
back  into  the  polar  cap.  Sojka  et  al.  (1986)  developed 
a  mathematical  model  with  a  Kr  dependence  that  can 
be  adjusted  to  produce  two.  three  or  four  cells  depend¬ 
ing  on  the  IMF  conditions.  The  above  models  con¬ 
centrated  on  explaining  phenomena  on  the  dayside  of 
the  polar  cap.  and  these  patterns  were  loosely  closed 
or  left  open  on  the  nightside.  Kan  and  Burke  (1985) 
tried  to  address  the  nightside  closure  while  still  keep¬ 
ing  the  four-cell  pattern  by  invoking  instabilities  and 
wa\c-likc  structures  in  the  nightside.  However,  none 
of  the  above  three-  and  four-cell  patterns  easily 
address  the  persistence  of  the  Harang  discontinuity 
during  both  quiet  and  B.  North  conditions  (Maynard. 
1974:  Heppnerand  Maynard,  1987). 

Hcppncr  and  Maynard  (1987)  empirically  derived 
typical  patterns  for  northward  IMF  conditions  for 
selected  passes  from  the  DE-2  electric  field  data  set 
(Maynard  et  al.,  1981)  in  which  the  IMF,  as  deter¬ 
mined  from  1 SEE  3.  was  clearly  unambiguous.  These 
patterns  represented  distortions  of  the  two-cell 
patterns  derived  for  southward  IMF  conditions 
(Hcppncr  and  Maynard,  1987).  The  distortions  were 
such  that  the  evening  convection  cell  rotated  into  or 
around  the  morning  cell.  The  distortion  increased  with 
an  increase  in  the  magnitude  of  the  northward  IMF 
component.  In  addressing  the  question  of  how  the 
convection  patterns  close  on  the  nightside.  the  dis¬ 
torted  two-cell  pattern  was  the  only  one  that  consis¬ 
tently  matched  the  data  for  the  large  majority  of  cases. 

This  poses  a  dilemma  in  deciding  which  convection 
pattern  is  appropriate  when  based  on  any  one  orbit 
of  data.  For  instance.  Killeen  et  al.  (1985),  following 
the  Burke  et  al.  (1979)  and  Potemra  et  al.  (1984) 
interpretation  of  a  “W”  pattern  as  indicative  of  four 
ceils,  explained  a  scries  of  passes  during  strong  B. 
North  conditions  in  terms  of  four-cell  patterns. 
Heppner  and  Maynard  ( 1 987)  showed  that  those  same 
passes  also  fit  the  distorted  two-cell  models.  It  would 
be  desirable  to  be  able  to  have  additional  tests  to  aid 
in  the  selection  of  the  best  interpretation. 

Ihe  Utah  State  Ionospheric  Model  (see  Schunk, 
1988)  provides  a  global  description  of  the  ionosphere 
based  on  input  parameters,  one  of  which  is  a  global 
ion  convection  pattern.  The  model  plasma  densities 
are  sensitive  to  the  past  convective  history  of  a  flux 
tube.  A  flux  tube  which  does  not  encounter  a  source 
of  ionization  such  as  the  sunlit  dayside  or  the  auroral 
region  will  quickly  become  depleted.  Thus,  convection 
patterns  which  have  central  portions  of  cells  which 
only  circulate  within  the  polar  cap  will  cause  electron 


density  holes  to  develop  in  the  polar  cap.  Evidence  of 
these  holes  will  start  to  develop  within  0.5  h  of  a 
pattern  change.  Sojka  and  Schunk  ( 1 987).  in  a  study  of 
the  response  of  the  ionosphere  to  B.  North  convection 
patterns,  caused  two  electron  density  holes  to  develop 
by  driving  the  ionospheric  model  with  a  theoretical 
four-cell  pattern.  By  applying  this  approach  using 
different  convection  patterns  as  a  driver,  one  can 
obtain  predictions  of  the  election  density  patterns 
in  the  high  latitude  ionosphere.  These  density  pre¬ 
dictions  can  then  be  qualitatively  compared  with 
measurements  of  electron  density  such  as  those  from 
the  DE-2  Langmuir  probe.  In  this  way  a  "blind"  test 
of  the  applicability  of  a  particular  convection  pattern 
can  be  performed  with  the  density  data  being  used  as 
the  discriminator. 

The  purpose  of  this  paper  is  (I)  to  test  whether 
distorted  two-cell  patterns  or  multiecl!  patterns  are 
most  applicable  under  relatively  stable  B.  North  input 
conditions  using  the  Utah  Stale  ionospheric  model 
and  DE-2  Langmuir  probe  data,  and  (2)  to  demon¬ 
strate  that  such  a  blind  test  can  be  used  to  restrict  the 
options  for  convection  pattern  configurations  for  a 
given  situation.  For  this  test  the  Utah  State  iono¬ 
spheric  mode!  was  driven  with  the  Heppner-Maynard 
(1987)  distorted  two-cell  model  and  with  the  Sojka 
et  al.  (1986)  three-cell  model  as  an  example  of  multi- 
cells.  No  other  input  conditions  were  varied.  The 
Sojka  et  al.  model  was  chosen  to  keep  the  similarity 
in  the  patterns  as  close  as  possible  as  explained 
below.  Times  were  chosen  for  comparison  with  the 
simulations  in  which  the  IMF  was  stably  configured 
and  when  DE-2  data  were  available.  Forty-one 
cases  were  considered.  The  results  of  this  test, 
which  favor  the  distorted  two-cell  model,  will  be 
generalized  in  the  Discussion  to  conclude  that  for 
individual  satellite  passes  the  prcscncc/abscncc 
and  position  of  polar  electron  density  holes  can  be 
used  to  discriminate  against,  or  rule  out.  specific 
convection  patterns.  As  an  example,  the  test  is 
applied  to  the  patterns  presented  in  the  theta  aurora 
paper  of  Frank  et  al.  (1986). 

IONOSPHERIC  MODF.I. 

The  ionospheric  model  adopted  for  this  study  was 
originally  developed  as  a  mid-latitude  multi-ion 
(NO+.  O^,  NT.  O*)  model  by  Schunk  and  Walker 
( 1 973),  and  was  later  extended  to  include  high  latitude 
cff.cts  due  to  convection  electric  fields  and  particle 
precipitation  by  Schunk  et  al.  ( 1975).  A  further  exten¬ 
sion  of  the  model  to  include  the  minor  ions  N*  and 
He+,  an  updated  photochemical  scheme,  and  the 
MSIS  atmospheric  model  (Hedin  et  al..  1977)  is 
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described  in  Schunk  and  Raitt  ( 1980).  More  recently, 
the  ionospheric  model  has  been  extended  by  Schunk 
and  Sojka  (1982)  and  Schunk  et  al.  (1986)  to  include 
thermal  conduction  and  diffusion-thermal  heal  How. 
so  that  the  ion  and  electron  temperatures  are  now 
rigorously  calculated  at  all  altitudes  between  120  and 
800  km.  ’ 

In  the  high  latitude  model,  tlux  tubes  of  plasma  are 
followed  as  they  convcct  through  a  moving  neutral 
atmosphere.  For  each  flux  lube,  altitude  profiles  of 
the  ion  and  electron  temperatures  and  the  NO".  0:~. 
N:*.  N".  O'",  and  He"  densities  arc  obtained  by  solv¬ 
ing  the  appropriate  continuity,  momentum,  and 
energy  equations  including  all  of  the  high  latitude 
processes  thought  to  be  important.  The  equations  are 
solved  over  the  altitude  range  from  120  to  800  km, 
with  boundary  conditions  specified  at  the  lower  and 
upper  ends  of  the  flux  tubes.  For  densities,  chemical 
equilibrium  is  assumed  at  120  km  and  no  escape  flux 
is  assumed  at  800  km.  For  temperatures,  local  thermal 
coupling  is  assumed  at  120  km  and  a  specified  heat 
flux  is  assumed  at  800  km.  For  this  study,  the  electron 
and  ion  heat  fluxes  through  the  upper  boundary  were 
taken  to  be  zero. 

Three-dimensional  density  and  temperature  distri¬ 
butions  are  obtained  by  following  many  flux  tubes  of 
plasma  and  by  keeping  track  of  their  locations  as  a 
function  of  time.  In  this  study,  all  runs  exceeded  24  h 
in  duration  to  ensure  that  equilibrium  was  reached 
for  the  normal  diurnal  variation  for  each  of  the  cases 
considered. 

The  ionospheric  model  does  not  self-consistently 
couple  to  the  magnetosphere  and  thermosphere  and. 
hence,  a  number  of  inputs  are  required  for  a  given  run. 
The  model  requires  global  distributions  for  the  neutral 
densities,  temperature,  and  wind.  At  high  latitudes, 
the  model  also  requires  the  magnetospheric  (con¬ 
vection)  electric  field  pattern  and  a  description  of  the 
auroral  electron  energy  flux.  These  important  thermo¬ 
spheric  and  magnetospheric  inputs  vary  from  one 
study  to  another:  those  adopted  for  this  study  arc 
described  in  the  following  two  sections  of  the  paper. 
A  complete  description  of  the  USU  Ionospheric 
Model  .s  given  by  Schunk  (1988)  and  the  reader  is 


referred  to  this  article  for  further  details.  A  com 


prehensive  review  of  the  model-data  comparisons  and 
model  applications  that  have  been  conducted  to  date 


is  given  by  Sojka  (1989). 


IONOSPHERIC  INPUTS 

Since  this  study  is  designed  to  contrast  two  types  of 
convection  models  and  their  effects  on  the  ionosphere. 


the  other  input  parameters  in  the  model  were  set  for 
specific  solar  cycle,  seasonal  and  magnetic  conditions 
and  held  constant  for  the  runs  with  each  convection 
pattern.  Table  1  lists  these  parameters.  For  the  neutral 
atmosphere  the  MSIS  model  (Hedin  et  ul..  1977)  was 
used,  while  for  the  neutral  winds  a  simple  pattern  with 
anti-sunward  flow  across  the  polar  cap  with  a  speed 
of  about  200  m  s  '  was  adopted.  This  pattefn  has  been 
used  in  several  previous  studies  with  the  ionospheric 
model  and  is  described  in  Sojka  et  al.  ( 1981). 

Electron  precipitation  in  the  auroral  oval  acts  as  a 
plasma  production  source,  a  source  of  bulk  electron 
heat,  and  a  source  of  energy  that  flows  through  the 
model’s  upper  boundary.  The  empirical  precipitation 
model  developed  by  Hardy  et  al.  (1985)  has  been  used. 
The  magnetic  activity  (Kr)  used  as  an  input  in  the 
Hardy  et  al.  model  is  given  in  Table  1 .  The  convection 
and  precipitation  patterns  were  chosen  such  that  they 
overlapped  as  closely  as  possible.  The  procedure  for 
calculating  the  ionization  rate  and  the  thermal  elec¬ 
tron  bulk  heating  rate  from  the  auroral  electron 
energy  flux  is  described  by  Schunk  (1988). 

CONVECTION  MODELS 

In  order  to  make  the  test  as  "blind”  as  possible,  the 
convection  models  selected  were  not  keyed  to  actual 
data,  but  were  based  on  the  IMF  conditions  at  the 
time  of  a  DE-2  pass.  In  the  high-latitude  ionosphere, 
the  convective  history  of  a  flux  tube  is  a  major  influ¬ 
ence  for  determining  its  electron  density.  Thus,  the 
total  convection  pattern,  not  the  detaiL  along  the 
satellite  pass,  becomes  the  controlling  factor.  Figure 
1.  for  IMF  northward  and  5,  positive,  shows 
the  electric  potentials  associated  with  the  Northern 
Hemisphere  distorted  two-cell  model  of  Heppner  and 
Maynard  (1987),  left  panel,  and  that  of  the  three-cell 
model  of  Sojka  et  al.  (1986).  right  panel,  which  is 
chosen  as  a  representative  of  multiple  cell  models.  The 
potential  is  represented  by  contours  drawn  at  2  kV 
intervals  in  a  magnetic  latitude-M.L.T.  coordinate 
system.  Negative  potentials  are  shown  as  dashed 
curves.  Figure  2  shows  the  corresponding  Northern 
Hemisphere  patterns  for  IMF  northward  and  B ,  nega- 
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auroral  regions  coincide.  The  three-cell  model  was 
chosen  to  provide  a  path  of  strong  anti-sunward  con¬ 
vection  on  the  appropriate  side  of  the  polar  cap  in 
order  to  conform  most  closely  with  its  accompanying 
distorted  two-cell  model.  The  rest  of  the  polar  cap 
was  allowed  to  vary  according  to  the  model. 

Differences  between  the  two  convection  models  are 
evident,  especially  in  the  noon  and  midnight  sectors ; 
however,  certain  similarities  exist.  Strong  anti-sun- 
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Tabll  I.  USU  lososi'M ric  Moon.  inputs 


Inpul  parameter 

Winter 

Summer 

Year  month/dav 

22  Dec.  1982 

22  June  1982 

Solar  1-10.7 

210 

210 

Magnetic  activity  Kr 

1 

1 

Magnetic  activity  Ar 

8 

8 

Auroral  precipitation 

Hardy  oval  (Kr  -  1) 

Hardy  oval  (Kr  =  1) 

Electric  field  distorted  2  cells 

Heppner  and  Mavnard 

Heppner  and  Mavnard 

Electric  field  mullicells 

Sojka  el  al. 

Sojka  et  al. 

Northward  IMF  B. 

15 

15 

B, 

±15 

Z  15 

ward  convection  occurs  at  the  dusk  side  of  the  polar 
e«.,i  for  B,  negative  and  at  the  dawn  side  for  B,  positive 
and  the  auroral  zones  are  similar.  Note  that  while  the 
Sojka  et  al.  (1986)  model  has  a  convection  cell  that 
circulates  only  entirely  within  the  polar  cap.  the  dis¬ 
torted  two-cell  equipotentials  return  to  the  auroral 
oval  for  each  circulation.  This  also  results  in  a  stronger 
sunwai  J  How  within  the  polar  cap  in  the  three-cell 
models  than  in  the  distorted  two-cell  models.  The 
selection  of  the  Sojka  et  al.  three-cell  model  as  the 
representative  of  multicell  models  was  done  to  max¬ 
imize  the  similarities  between  the  distorted  two-cell 
models  and  the  multicell  models.  More  cells  within 
the  polar  cap  would  have  increased  the  differences 
between  the  multicell  and  the  Heppner-Maynard  dis¬ 
torted  two-cell  models  and.  as  shown  by  Sojka  and 
Schunk  (1987).  would  have  increased  the  number  of 
polar  holes  created  by  the  model  in  the  electron 
density  patterns. 

RESUI.TS  ANI)  COMPARISONS 

A  total  of  41  satellite  crossings  of  the  polar  region 
which  had  stable  IMF  conditions,  were  selected  for 
study.  The  DE- 2  passes  were  selected  from  those  used 
by  Heppner  and  Maynard  (1987)  in  the  derivation  of 
the  distorted  two-celi  models  and  represent  cases  in 
which  the  IMF  direction  as  determined  from  /SEE  3 
data  was  stably  configured  from  I  5  h  before  the  pass 
to  30  mm  before  the  pass  to  ensure  that  the  IMF 
conditions  wuc  unambiguous.  These  passes  lep- 
resented  a  readily  accessible  data  set.  Any  individual 
pass  can  be  represented  by  many  solutions  for  the 
global  potential  pattern.  Because  the  test  did  not  fine 
tune  the  models  to  individual  passes  and  because  the 
Heppner-Maynard  models  are  meant  to  be  average 
patterns  and  not  represent  any  one  pass  exactly,  the 
use  of  these  passes  does  not  introduce  a  significant 
bias  in  the  results.  The  initial  scan  used  appropriate 
winter,  summer  or  equinox  conditions  in  which  the 


ionospheric  model  was  run  to  equilibrium.  The  appro¬ 
priate  universal  time  was  then  selected  for  the  model- 
to-data  comparisons.  For  a  selected  group  of  orbits 
the  model  was  rerun  to  match  more  closely  the  sun¬ 
spot  number  and  convection  pattern  size.  Examples 
presented  in  the  paper  are  from  that  subset. 

Before  presenting  model-data  comparisons,  it  is 
important  to  note  that  the  ionospheric  model  was  run 
to  obtain  diurnally  reproducible  density  features.  In 
this  case,  the  ionosphere  has  sufficient  time  to  fully 
adjust  to  the  driving  convection  pattern  and.  hence, 
the  distinguishing  density  features  associated  with  the 
convection  pattern  will  be  well  developed.  In  practice, 
however.  B.  will  not  remain  steady  for  a  24-n  period. 
Nevertheless,  when  the  IMF  changes  direction,  the 
time  constant  for  corresponding  density  changes  is  of 
the  order  of  0.5  h.  It  is  within  this  time  period  that  the 
distinguishing  density  features  should  start  to  appear, 
although  they  may  not  be  as  well  developed  as  the 
equilibrium  solutions  predict.  For  this  reason,  and 
also  considering  that  the  convection  patterns  are  not 
exact  in  the  representation  of  actual  conditions,  only 
a  qualitative  comparison  between  model  predictions 
and  measurements  is  possible. 

Figure  3  presents  a  comparison  of  results  from  a 
winter  case  in  which  B,  was  positive.  The  model 
densities  for  the  different  convection  patterns  are  pre¬ 
sented  on  the  left.  The  top  polar  dial,  encompassing 
from  50  magnetic  latitude  to  the  pole,  displays  oxy¬ 
gen  ion  densities  at  800  km  (approximately  equivalent 


to  total  density)  for  the  Ilcppncr  Maynard  pattern 
run.  while  the  bottom  polar  dial  shows  the  oxygen  ion 
density  for  the  muiticell  pattern  run.  The  DE- 2  track 
is  given  by  the  dark  line  across  the  plot  from  dawn  to 
dusk.  The  measured  electron  density.  Nc ,  and  tem¬ 
perature.  Tt  (Krehbicl  et  al..  1981),  and  horizontal 
electric  field  along  the  orbit  track.  £,  (Maynard  et 
al.,  1981),  from  DE- 2  is  shown  on  the  right.  The 
altitude  for  this  pass  is  near  700  km  in  the  center  of  the 
polar  cap.  The  comparisons  were  done  qualitatively 
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Electric  potential  contours  for  the  Northern  Hemisphere  associated  with  the  Heppner  and  Maynard  (1987),  left  panel,  and  the  Sojka 
el  al.  (1986),  right  panel,  models  for  IMF  northward  and  B,  positive  conditions. 

The  potentials  are  contoured  at  2  kV  intervals  in  a  magnetic  latitude-M.L.T.  coordinate  system. 
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between  the  global  patterns  and  the  density  traces. 
However,  in  order  to  show  clearly  the  differences  and 
similarities,  the  heavy  solid  and  dotted  lines  in  the 
density  plot  arc  actual  values  along  the  orbit  track 
calculated  from  the  ionospheric  model  runs  using  the 
Heppncr-Maynard  distorted  two-cell  and  multicell 
models,  respectively.  The  model  results  using  the  multi- 
cell  pattern  show  a  deep  hole  in  the  center  of  the 
polar  cap.  Both  model  runs  have  an  enhanced  morn¬ 
ing  side;  however  the  Heppner-Maynard  distorted 
two-cell  tun  does  not  have  a  deep  hole  in  the  polar 
cap.  showing  only  a  small  reduction  at  the  edge  of  the 
auroral  zone.  Correspondingly,  the  measured  electron 
density  is  highest  on  the  morning  side  of  the  polar 
cap.  with  the  high  density  region  extending  to  the 
noon-midnight  meridian.  This  is  the  region  of  the 
highest  anti-sunward  convection  speeds  both  in  the 
data  and  in  the  models,  bringing  highly  populated  flux 
tubes  from  the  sunlit  daysidc  region  to  the  polar  cap. 
The  lack  of  a  significant  hole  in  the  measured  density 
at  the  pole  is  obvious.  In  this  respect  the  measured 
density  profile  best  fits  the  density  results  obtained 
using  the  Heppner-Maynard  convection  model. 

Figure  4  shows  an  example  in  the  same  format  for 
fl,  negative.  In  this  case  the  satellite  altitude  at  the 
center  of  the  polar  cap  is  near  400  km.  The  polar  dials 
depict  the  electron  density  at  300  km  from  the  model 
runs.  Note  ihat  both  convection  patterns  produce 
polar  holes.  For  the  distorted  two-cell  pattern,  the 
hole  is  centered  on  the  morning  side  of  the  noon- 
midnight  meridian,  while  for  the  multiccll  pattern  the 
hole  is  centered  slightly  to  the  evening  side  of 
the  noon-midnight  meridian.  The  measured  electron 
density  profile  in  Fig.  4  shows  a  deep  hole  on  the 
morning  side  of  the  noon-midnight  meridian  plus  a 
smaller  hole  on  the  equatorward  edge  of  the  evening 
auroral  region  (seen  in  both  models).  This  lower  lati¬ 
tude  hole  is  the  midlatitude  trough  (Spiro  el  ul..  1978). 
The  electric  field  in  this  case  is  more  irregular,  but  still 
has  an  overall  pattern  consistent  with  the  distorted 
two-cell  pattern.  The  density  data  fit  the  distorted  two¬ 
cell  driven  results  better,  as  can  be  seen  by  comparison 
with  the  heavy  lines  which  show  the  model  results 
along  the  orbit  track. 

The  two  examples  favor  the  Heppner-Maynard  dis¬ 
torted  two-cell  patterns  very  strongly.  Each  of  the  4! 
passes  was  compared  with  the  model  results,  ar.d  the 
agreement  was  categorized  as  being  strongly  in  favor 
of.  moderately  in  favor  of.  or  weakly  in  favor  of  the 
appropriate  convection  driver.  Two  other  categories 
were  required.  First,  several  passes  were  fitted  equally 
well  by  both  models  and  were  categorized  as  "no 
choice".  In  several  of  the  passes  for  Bs  negative, 
neither  convection  pattern  provided  reasonable  com¬ 


parisons  and  these  passes  were  judged  as  “other". 
Some  of  the  passes  categorized  as  “other"  were  passes 
where  the  size  of  the  auroral  oval  and  polar  cap  was 
much  larger  than  that  used  in  the  models  due  to  an 
extremely  large  total  magnitude  of  the  IMF.  In 
addition,  lack  of  full  agreement  in  many  cases  could 
be  expected,  since  only  the  most  distorted  convection 
patterns  were  used  in  the  ionosphere  model  runs  even 
though  the  less  distorted  patterns  might  in  some  cases 
have  been  more  applicable.  The  results  arc  presented 
in  a  bar  graph  in  Fig.  5.  A  full  spectrum  of  cases  is 
seen ;  however,  a  majority  of  the  cases  (23)  favor  to 
some  degree  the  Heppner-Maynard  distorted  two- 
ccll  model.  Ten  of  the  41  passes  considered  strongly 
favored  the  distorted  two-cell  patterns.  When  the  no 
decision  cases  arc  omitted,  the  results  are  even  more 
strikingly  in  favor  of  the  distorted  two-cell  model 
(23  of  29). 

DISCUSSION 

To  put  the  above  results  in  perspective,  this  test 
only  compared  the  Heppner-Maynard  distorted  two¬ 
cell  model  against  the  Sojka  el  al.  three-cell  model. 
The  passes  studied  were  from  the  group  which  had 
stable  IMF  conditions  and  had  been  used  to  derive 
the  distorted  two-cell  patterns  (as  explained  earlier, 
this  was  not  considered  a  biasing  factor).  While  the 
results  were  clearly  in  favor  of  the  distorted  two-cell 
model  used,  they  do  not  a  posteriori  favor  all  distorted 
two-cell  models  over  all  multicell  models.  However, 
certain  characteristics  of  the  comparisons  do  lend 
themselves  to  generalization. 

The  primary  distinguishing  -haracteristics  used  to 
categorize  the  agreement,  or  lack  thereof,  were  the 
depth  and  position  of  polar  cap  holes  in  the  electron 
density.  The  polar  cap  holes  arc  produced  in  the  iono¬ 
spheric  model  when  magnetic  flux  tubes  circulate  only 
within  the  polar  cap  (i.e.  on  equipotential  contours 
that  are  confined  to  the  polar  cap),  and  particularly 
when  these  flux  tubes  remain  in  darkness.  In  the  -B, 
case  the  distorted  t<  -cell  pattern  develops  a  hole 
because  of  the  weak  ionization  levels  encountered 
when  the  circulating  flux  lubes  contact  the  post-mid¬ 
night  auroral  oval.  The  holes  are  deepest  at  universal 
times  when  the  terminator  is  farthest  toward  the  day¬ 
sidc  of  the  magnetic  pole.  The  choice  between  the 
convection  patterns  depends  on  the  location  of  the 
hole  and  is.  therefore,  model  sensitive.  For  instance, 
the  Sojka  and  Schunk  ( 1987)  four-cel!  pattern  created 
two  holes  in  the  center  of  the  polar  cap  because  two 
of  the  cells  had  centers  that  only  circulated  in  the  polar 
cap.  Different  distorted  two-cell  or  multicell  patterns 
could  result  in  changes  in  the  location,  size  and  or 
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Fig  5.  Statistical  comparison  of  DE-2  electron  density  data  with  model  results  for  a  northward 

IMF  WITH  B,  POSITIVE  (TOP  PAN'LL)  AND  B,  NEGATIVE  (BOTTOM  PANEL). 

The  seven  histograms  associated  with  each  pane!  represent,  from  left  to  right,  the  following  conditions . 
strongly  in  favor  of.  moderately  in  favor  of.  weakly  in  favor  of  the  Hcppncr-Maynard  (1987)  model, 
equally  likely:  and  then  weakly  in  favor  of.  moderately  in  favor  of.  and  strongly  in  favor  of  the  Sojka  ct 
al.  (1986)  model.  Of  a  total  of  41  data  sets,  six  could  not  be  grouped  according  to  these  descriptions  and. 

hence,  arc  in  the  "other"  group. 


20  PASSES 


existence  of  holes  and.  therefore,  changes  in  the  com¬ 
parison  statistics. 

Thus,  the  presence  or  absence  and  the  position  of 
holes  in  the  electron  density  in  the  polar  cap  can  be 
used  as  a  criterion  to  rule  out  possible  convection 
patterns.  For  example,  in  Fig.  3  for  Northern  Hemi¬ 
sphere  B,  ■ .  the  distorted  two-cell  pattern  docs  not 
produce  a  clear  hole,  whereas  the  multicell  pattern 
does.  No  hole  is  observed.  Any  other  convection 
models  that  have  large  cells  whose  centers  only  cir¬ 
culate  within  the  polar  cap  will  not  fit  the  observed 
data.  Also  the  dawn  side  electron  densities  in  the  polar 
cap  require  significant  anti-sunward  flow  from  the 
highly  ionized  dayside  m  this  B ,  -  case.  Both  of  these 
features  are  more  easily  satisfied  by  the  distorted  two¬ 
cell  patterns. 

The  case  that  most  favored  the  multicell  pattern  in 
the  B ,  data  set  (not  shown)  involved  a  condition  in 
which  B,  was  nearly  zero  and  the  principal  horizontal 
component  was  in  Bx.  The  distorted  two-cell  pattern 
breaks  down  in  this  situation  of  B.  »  either  from 
the  distortions  during  transition  between  the  patterns 
or  from  a  different  forcing  function.  This  IMF  con¬ 
dition  produced  a  four-cell  pattern  in  the  study  of 
Clauer  and  Friis-Christensen  (1988). 


be  adapted.  These  were  generally  cases  where  the 
patterns  were  not  sized  correctly  due  to  an  abnormally 
high  B.  or  cases  where  the  potential  data  were 
extremely  variable  and  bore  no  resemblance  to  the 
model  patterns.  Thus,  no  one  set  of  patterns  will  rep¬ 
resent  all  conditions.  However,  the  distorted  two-cell 
concept  offers  a  new  approach  to  changing  from  B. 
South  to  B.  North  conditions  in  a  smooth  manner.  It 
does  not  cover  the  case  when  the  electric  field  is  so 
variable  that  no  background  convection  pattern  can 
be  discerned.  Hoffman  el  al.  (1988)  found  for  quiei 
conditions  with  B.  North  that  a  pattern  was 
recognizable  in  the  electric  field  60%  of  the  time,  even 
in  their  generally  weakly-driven  cases.  Also,  it  should 
be  noted  that  during  northward  IMF  conditions  some 
precipitation  is  present  in  the  polar  cap.  However,  for 
the  ionospheric  model  runs,  polar  cap  precipitation 
was  not  included  because  it  was  not  clear  how  to 
incorporate  it  in  the  Hardy  et  al.  (1985)  precipitation 
model. 

The  above  generalization  restricting  the  topology 
of  convection  patterns  based  on  the  existence  and 
position  of  holes  in  the  electron  density  in  the  polar 
cap  can  be  applied  to  the  theta  aurora  passes  presented 
by  Frank  et  al.  (1986)  to  test  the  multicell  patterns 
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Fig.  5.  DF.-2  electron  density  plots  for  three  or  the  theta  aurora  passls  discussed  by  Frank  <>/ «/ 

(19X61 

The  black  bar  on  the  horizontal  axis  of  each  plot  denotes  the  location  of  the  auroral  emission  identified 

as  the  theta  bar. 
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passes  were  identified  as  having  multicell  patterns. 
Passes  1437  (8  November  1981)  and  1323  (31  October 
1981)  were  for  B,.  negative  conditions.  Note  that  for 
pass  1437  9.  was  only  weakly  positive.  For  pass  1111 
(17  October  1981)  the  B,  component  was  variable. 
Figure  6  shows  the  LZ-2  electron  density  plots  for 
these  passes.  Only  orbit  1 323  shows  a  clear  polar  hole. 
The  hole  is  somewhat  to  the  morning  side  of  the 
noon-midnight  line.  In  comparing  with  Fig.  4  and  the 
discussion  of  that  figure,  this  would  favor  the  distorted 
two-cell  pattern  over  the  three-cell  pattern.  Frank  et 
al.  (1986:  Figs  II  and  7)  show  a  polar  cell  on  the 
evening  side  and  an  indeterminate  structure  on  the 
morning  side.  Since  their  evening  cell  could  not  pro¬ 
duce  the  observed  hole,  this  test  rules  out  their  pattern . 
Pass  1 1 1 1  has  almost  no  discernable  density  structure. 
It  tracks  from  about  9  to  21  h,  passing  through  holes 
for  either  B,  condition  for  the  multicell  runs  in  Figs  3 
and  4.  The  orbit  track  would  just  touch  the  top  of  the 
hole  for  B,.  negative,  and  there  is  no  hole  for  B,  positive 
in  the  distorted  two-cell  runs.  The  lack  of  a  deep  hole 
could  be  attributed  to  the  near  equinox  conditions 
which  would  result  in  more  of  the  cap  being  in 
sunlight.  However,  even  in  summer,  the  multicell  con¬ 
vection  driven  model  run  showed  some  depletion.  The 
lack  of  density  structure  is  inconsistent  with  the  multi¬ 
cell  pattern  proposed  by  Frank  et  al.  (1986:  Figs  12 
and  9).  For  pass  1437  no  large  polar  hole  is  seen; 
however,  dips  of  the  order  of  a  factor  of  two  or  less 
are  seen  in  the  evening  side  of  the  polar  cap  and  in 
the  vicinity  of  the  theta  bar.  The  Sojka  and  Schunk 
(1987)  four-cell  hypothetical  model  for  a  theta  aurora 
condition  had  a  ridge  (not  a  dip)  separating  the  two 
polar  holes  partly  due  to  particle  production  from 
electron  precipitation  at  a  simulated  theta  aurora 
location.  This  pass  does  not  agree  well  with  cither  of 
the  simulations  in  this  test,  although  there  is  certainly 
no  evidence  of  having  crossed  cells  with  circulation 
only  within  the  polar  cap,  as  could  be  inferred  from 
the  pattern  for  this  pass  given  by  Frank  et  al.  (1986: 
Figs  10  and  5).  The  fact  that  one  of  the  small  dips 
occurs  at  the  region  of  the  precipitation  is  contrary  to 
both  the  Frank  et  al.  and  the  Sojka  and  Schunk 
patterns.  The  density  configurations  from  these  three 
passes  suggest  that  the  multicell  patterns  proposed  by 
Frank  et  al.,  which  would  have  produced  polar  holes 
on  the  evening  side  of  the  noon-midnight  meridian, 
were  not  correct  for  those  conditions.  In  all  three  of 
the  above  passes  this  test  does  not,  however,  define 
the  correct  convection  configuration. 

CONCLUSIONS 

This  “blind"  test  of  convection  patterns  using  the 
ionospheric  mode!  demonstrates  a  preference  for  dis¬ 


torted  two-cell  patterns  in  the  majority  of  cases 
studied.  Preference  was  determined  by  comparing  the 
presence  (or  absence)  and  location  of  polar  holes  in 
measured  data  vs  model  predictions.  While  this  test 
does  not  uniquely  specify  the  best  pattern,  preference 
for  distorted  two-cell  patterns  rejects  convection  pat¬ 
terns  which  only  circulate  within  the  polar  cap  for 
B,.  conditions  and  restricts  shapes  for  B,.-  conditions. 
As  a  corollary,  convection  patterns  whose  center  cir¬ 
culates  only  within  a  dark  polar  cap  must  produce 
holes  in  the  electron  density.  This  means  that  electron 
density  profiles  can  be  used  to  exclude  certain  ad  hoc 
convection  models  proposed  to  explain  the  data  from 
a  single  pass ;  however,  they  cannot  be  used  to  derive 
unambiguously  the  correct  model.  Better  statistics  arc 
needed  for  future  evolution  of  patterns  for  B.  North, 
whether  they  be  of  the  distorted  two-cell  or  multiccll 
variety, 
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